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The atomic parity non-conservation(PNC) experiments has reached accuracies which have im- 
portant implications for physics beyond the standard model. An optical rotation experiment to 
measure the atomic PNC of the 3 Pi(6s6p) — 1 Pi(6s6p) in Yb was proposed recentlyjlj. Our screened 
electron-electron coulomb potential multi-configuration Dirac-Fock calculation of th PNC induced 
El transition amplitude of this transition EIpnc = —96.0 x 10~ n ieao(— Qw/N) of 171 Yb is more 
than two orders of mag nitude larger than £l PNC (6s-7s) = -0.8991 x 10' 11 iea (-Q w /N) of 133 Cs. 



I. INTRODUCTION 

The nuclear weak charge Qw of an atom is a parameter 
of the atomic PNC arising from the nucleon-electron vec- 
tor axial vector interaction. It is possible to extract Qw 
by measuring the PNC induced electric dipole transition 
amplitude -EIpnc and combining it with atomic theory 
calculations. The .EIpnc of the 133 Cs 6s — 7s transi- 
tion has been measured to very high accuracy Q. The 
Qw of 133 Cs obtained from the experimental data after 
combining with the atomic theory calculations || ^| was 
further refined to 72.06(28) cxp t(34)thoor after the mea- 
surement of the 6s — 7s transition polarizabilityQ. This 
has a 2.5<7 deviation from the predictions of the stan- 
dard model Qw = 73.20(13) ||. The uncertainty of the 
experiment results and atomic theory calculations are es- 
timated at .35% and 1% respectively. This indicate the 
possibility of improving the results further by reducing 
the atomic theory uncertainty. 

The recent atomic theory calculations which treat the 
Breit interaction more rigorously have introduced correc- 
tions of .9%|), .6%jg and .4%g to Q w of 133 Cs. The 
variation of the results is due to the difference of the cal- 
culation methods and many-body effects included. The 
uncertainties of these calculations are estimated at 1%. 
Considering the unique implications of the parameters 
obtained from the atomic PNC phenomena to the physics 
beyond the standard model, it is desirable to reduce the 
atomic theory uncertainties and also confirm the 133 Cs 
results using other atoms. 

An important criteria of choosing a candidate atom is 
the Z 3 a scaling of the PNC interaction, which indicates 
the advatage of choosing high Z atom. The 1 So — 3 D\ 
transition of atomic Yh(Z = 70), which has been stud- 
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ied theoretically 0, |TT|, [[| || and experimentally |l5| 
|27[ in detail is a promising candidate of the ongo- 
ing and future atomic PNC experiments. This tran- 
sition is a suitable choice due to the nearly degener- 
ate 3 L»i(5d6s)(24489cm- 1 ) and 1 P 1 (6s6p)(25068cm~ 1 ), 
which can enhance the PNC mixing. The PNC mixing 
enhancement between this pair of levels can also be ex- 
ploited in the 3 Po — 1 Pi transition. 

The 3 Po — 1 Pi transition was proposed for atomic PNC 
experiment in a recent workjl) and estimates £Tpnc ~ 
7 x 10~ 10 eao which is almost one order of magnitude 
larger than that of 133 Cs -E1 PNC ~ 8 x 10 _11 ea . The use 
of the metastable state 3 Po as the initial state offers the 
possibility of separating only the atomic density depen- 
dent plane of polarization rotation angle from the wave- 
length dependent systematic errors. This is expected to 
better the results obtained from the Stark-interference 
experiments, where the use of the high intensity laser 
fields limits the statistical precision due to the light shift. 

In this paper we have investigated the £J1pnc of the 
3 Po.i — 1 P\ transition using multi-configuration Dirac- 
Fock method. The atomic theory calculations are nec- 
essary to extract the nuclear weak charge Qw from the 
experimental value of £Tpnc- To check the accuracy 
of the atomic calculations we also study the transition 
properties and hyper fine structure constants. All the 
calculations are in atomic units ( % = e = m e = 1). 



II. THE METHOD OF CALCULATION 

The nuclear spin-independent atomic PNC arises from 
the axial-vector vector elcctron-nucleon interaction com- 
ponent of the neutral weak current interaction between 
the electrons and nucleons of an atom. The interaction 
is mediated by Zq bosons, which is a prediction of the 
electron-weak unification. The effective form of the in- 
teraction Hamiltonian can be obtained by treating the 
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nuclear part non-relativistically 
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(2.1) 



where Gf is the Fermi coupling constant, 75 is the Dirac 
matrix and p nuc is the nuclear density. The PNC in- 
duced El transition amplitude between the initial and 
final states and f) is 
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(2.2) 



where |\E'/) are the intermediate states which are opposite 
in parity to and |^/), and Ei, Ef and Ei are the 
energies of the states. 



The numerator of (2.2) has Hp^ l c , effective only within 
the nuclear region due to the nuclear density p nuc (r) and 
the dipole D — —r in length gauge, which has significant 
contribution from the large radial range. And the ener- 
gies in the denominator has large contribution from the 
mid radial range where the electron density is high. The 
calculation of the I?lpfj C require atomic state functions 
which are accurate over all radial ranges. -EIpnc can be 
calculated once the atomic wave-functions are known. 



A. Wave-function Calculation 

A method suitable for rare earth atoms which has large 
configuration mixing is the multi-configuration Dirac- 
Fock (MCDF), which is the relativistic adaptation of the 
multi-configuration Hartree-Fock. The MCDF approxi- 
mates an atomic state function \TPJM) as a linear com- 
bination of configuration state functions (CSF) \jPJM), 
which are again a linear combination of Slater determi- 
nants. Where P, J and M are the parity, total angu- 
lar momentum and magnetic quantum numbers respec- 
tively, and r and 7 are the additional quantum numbers 
required to defined the ASF and CSFs uniquely. Then 



\TiPJM) =Y,c r r l \lrPJM), 



(2.3) 



where c r r ; are the coefficients of the CSFs. The en- 
ergy functional is defined using a set of ASFs which mix 
strongly and the Dirac-Coulomb Hamiltonian 
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where cti and (3i are the Dirac matrices, Pi is the momen- 
tum and N is the number of electrons. The orbitals are 
then generated variationally. The orbitals from the neg- 
ative continuum are also the solutions of such a method, 
however only the bound states can be chosen by impos- 
ing the boundary condition that the orbitals ip(r) — > as 



r — > 00. Choosing the ASFs contributing to the valence- 
valence correlation effects, the orbitals captures impor- 
tant correlation effects. 

Though MCDF method can represent the valence- 
valence correlation, it is not a suitable method to cal- 
culate core- valence and core-core correlation effects. A 
large number of CSFs is required to represent these cor- 
relation effects, which is computationally difficult as self 
consistent field method like MCDF is not suitable for a 
calculation involving large CSF spaces. However, these 
correlation effects can be calculated using configuration 
interaction(CI) using a set of virtual orbitals generated 
from the MCDF potential. 

A large set of CSFs is required in the CI calculation to 
represent all classes of correlation effects. A more efficient 
method is to model it using the most important CSFs 
and modifying the form of the electron-electron Coulomb 
interaction potential as 
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£^-§1^(1) -C^(2), 



(2.5) 



K 



where chk are the constants that modifies the K mul- 
tipole of the inter electron Coulomb interaction poten- 
tial. The suitable values of these constants can be ob- 
tained by matching properties like excitation energies 
calculated using this potential with the experimental 
data. The earlier calculations have shown that, this ap- 
proach can reproduce the experimental data to very good 
agreement [ flO[ [l6| . The values of a are chosen such that 
0. < oik < 1< and for consistency the screening parame- 
ters are used in the MCDF calculations to generate the 
orbitals. 



B. Properties Calculation 

The most important transition properties like oscilla- 
tor strength, life time, polarizability, etc arise from the 
electric dipole transition. And all of these depend on the 
reduce matrix elements of the electric dipole operator. 
Consider the initial and final states of the atom and 
then from Wigner-Eckert theorem 
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(2.6) 

where Jf and Ji are the total angular momenta of the fi- 
nal and initial states, m/ and are the magnetic quan- 
tum number of the final and initial states, and q is the 
component of the dipole operator. ( x f , /||_D||^E'i) = Dfi is 
the reduced matrix element and it is independent of ge- 
ometry. The reduced matrix elements can be expressed 
in terms of the ASFs as 

Dfi = (TfPfJfWDWTM) (2.7) 

= ^Cr^Csr^rPfJfWDlhsPiJi)- (2.8) 
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The oscillator strength of the transition 

2AE 



ffi 



3(2J/ + 1) 



fi\ 



\9 f ) is 
(2.9) 



where AE = Ef — Ei. The independent particle ap- 
proximation of the inter-electron electromagnetic inter- 
action implies that it is not gauge invariant. However, 
the gauge invariance is restored if the inter-electron elec- 
tromagnetic interaction is represented completely by in- 
cluding the correlation effects correctly. The agreement 
of the dipole matrix elements or transition properties cal- 
culated in different gauges indicate the completeness of 
the correlation effects included, which can also be inter- 
preted as the accuracy of the wave-functions. 

The hyperfine constants are the measure of the 
strength of the electron-nucleus parity allowed electro- 
magnetic multipole interactions. The general form of the 
interaction Hamiltonian is 
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where and M<- k "> are spherical tensor operators of 
rank k in the electron and nuclear space respectively. 
These represent the electromagnetic multipoles. The k = 
1 and 2 corresponding to the magnetic dipole and electric 
quadrupole respectively are the most important. The 
atomic states are then the eigenstates of the total angular 
momentum 



F = I + J, 



(2.11) 



where I is the nuclear spin. The shift in energy due 
to the magnetic dipole and electric quadrupole hyperfine 
interactions of the atomic state Iji'yjIJFMf) are 
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where C = F(F + 1) - J(J+ 1) - II- 1), and Aj and Bj 
are the magnetic dipole and electric quadrupole hyperfine 
constants respectively. These are defined as 
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where fj,i and Q are the nuclear magnetic dipole and elec- 
tric quadrupole moments. 



III. RESULTS 

A. Screened Coulomb Potential 

The low-lying levels of the Yb are given in the Fig 
[I]. The ifpNC mixing between the \ 1 Pi(6s6p)) and 



\ 3 Di(5d6s)) reduces to mixing between 5d and 6p at the 
single particle level, which is negligible as 5d is almost 
zero in the nuclear region. However, 6s and 6p mixing 
can occur through the mixing ( 3 Di(5d6s)|i/pNc|5d6p), 
which arises due to the strong configuration mixing be- 
tween |6s6p) and \5d6p). So the most important con- 
tribution to the -Epnc( 3 -Po — 1 -Pi) transition amplitude 
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( 1 P 1 (6s6p)\H PNC \5d6p)(5d6 P \D\ 3 P {6s6 P )} 



Ei 



5d5p) 



This indicates that the choice of | 3 Pi(6s6p)) as ini- 
tial state offers the possibility of non-zero contribution 
from the leading configurations to the normal and con- 
jugate terms. Whereas the largest contribution to the 
^pnc( 1 ^o(6s 2 )- 1 P 1 (6 S 6p)) is 



{ 1 P 1 (6s6p)\D\5d6p){5d6 P \H PNC \ 1 S (6s 2 )) 



Ei. 



S (6s6p)) 



Ei 
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(3.1) 



the conjugate term does not contribute as |5d6p) is dou- 
ble excitation with respect to j 1 Sq(6s 2 )) . 
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FIG. 1: The low- lying levels of atomic Yb. The double arrow 
represents the -Hpnc mixing and single arrow represents the 
transitions of interest for _E1pnc measurement. 

All the calculations are done using GRASP2[jl7j pack- 
age. The choice of the screening parameters as and the 
dependence of the excitation energies on these parame- 
ters are discussed in an earlier work|l(| [l(J. The values 
of the ax s which can give excitation energies very close 
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to the experimental data are ao = 0.997, a% = 0.667, 
« 2 = 0.980 and rest are set to unity. The even and odd 
parity CSF space consist of the following non-relativistic 
configurations : 
Even parity 

(5/4/ 14 ) (6s 2 + 5d6s + 5d 2 + 6p 2 )) 

+ (5/4/ 13 ) (6s 2 6p + 5d 2 6p + 5d6s6p) 

+ (5p 5 4/ 14 ) (6s 2 6p + 5d 2 6p + 5d6s6p) 

Odd parity 

(5/4/ 14 ) (6s6p + 5d6p) 

+ (5p 6 4f 13 ) (5d 2 6s + 5d6s 2 + 6s6p 2 + 5d6p 2 ) 

+ (5p 5 4/ 14 ) (5d 2 6s + 5cffis 2 + 6s6p 2 + 5d6p 2 ) 

The configuration space is chosen to represent the im- 
portant valence-valence and core-valence correlation ef- 
fects. The excitation energies calculated using these 
choices of a^s and configurations are given in Ta- 
ble!. ^ or a com P arat i ve study we have also calculated 



TABLE III: iJlpNC of recently studied and calculated 
atomic/ionic transitions. These are given in units of 
iea (-Qw/N) x 10" 11 



Atom/ion 


Z 


N 


Transition 
Initial state final state 


PIpnc 


Cs 


55 


78 


!6Sl/2> 


\7s 1/2 ) 


-0.8991(36)" 


Cs 


55 


78 


!6Sl/2> 


\5d 3/2 ) 


3.75 6 


Ba+ 


55 


82 


|6Sl/2> 


|5d 3/2 ) 


2.17 6 


Yb 


70 


101 


rSo(6s 2 )) 


j 3 £>i(5d6s)) 


79.38 c 


Tl 


81 


124 


|6pi/ 2 > 


|6p 3 /2> 


66.7± 1.7 d 


Fr 


87 


136 


|7Sl/2> 


!6d 3 / 2 ) 


57.1 b 


Ra 


88 


139 


rSo(7 S 2 )) 


j 3 £>i(7s6d)> 


77.0 e 


Ra 


88 


135 


r5o(7s 2 )) 


| 3 L>i(7s6d)) 


76.0 e 


Ra+ 


87 


135 


|7«l/2> 




42.9 b 




TABLE I: The excitation energies of the low-lying levels in 

cmT 1 



Level 


Expt. data 


Theory 


3 Pi(6s6p) 


17992 


17651 


3 P 2 (6s6p) 


19710 


19855 


3 Di(5d6s) 


24489 


24467 


1 Pi(6s6p) 


25068 


25024 



£1pnc( 3 -Pi(6s6p) - 1 Pi(6s6p)). However | 3 Pi(6s6p)) 
can decay to the ground state I 1 So (6s 2 )) by magnetic 
quadrupole transition. Where as | 3 Po(6s6p)) — | 1 S'o(6s 2 )) 
is highly forbidden. The calculated values of reduced ma- 
trix element of E'Ipnc are given in Table 0. For com- 



TABLE II: The PIpnc reduced matrix elements of atomic 
171 Yb(J = 1/2) in units of iea (-Q w /N) x 10 -11 . 



Transition 


| PlpNC | 


1 S (6s 2 )- 3 Di(5d6s) 


77.8 


3 Pi(6s6p) - 1 Pi(6s6p) 


-71.7 


3 P (6s6p) - 1 Pi(6s6p) 


-96.0 



parison ||_B1pnc( 1 'S'o(6s 2 ) — 3 £>i(5g?6s))| | is also included 
in the table. The ||£TpncII of atoms and ions which has 
been studied and calculated recently are given in Table 



III. The comparison of the iTpNc transition amplitudes 
indicates that the 3 Po(6s6p)— 1 Pi(6s6p) transition ofYb 
is the largest. 

B. MCDF Calculations 

The occupied orbitals (1 — 6)s, (2 — 5)p, (3 — i)d and 
4/, and the valence orbitals 6s, 5g? and 6p are gener- 



ated by a sequence of MCDF calculations and non of the 
orbitals are frozen in each of the calculations. The non- 
relativistic configurations of the final calculation of the 
sequence are: 6s 2 , 6p 2 , 5d6s, 5c? 2 , 6s6p and 5d6p. These 
orbitals are generated spectroscopic, that is the number 
of nodes satisfy the central field condition. This selection 
of non-relativistic configurations include the leading con- 
figurations of the 3 P!(6s6p), 3 £>i(5g?6s) and 1 P i (6s6p) 
levels, and the configurations which mix strongly with 
these. The levels obtained as a result of the MCDF cal- 



culation are given in the column MCDF1 of Tabic IV . The 



TABLE IV: The Yb levels in units of cm -1 calculated using 
MCDF method with the leading and next leading configura- 
tions 



Levels 




Excitation Energies(E) 






Expt 


MCDF1 


MCDF2 


3 Po 


17288 


13428 


13535 


3 Pi 


17992 


14119 


14234 


3 P 2 


19710 


15644 


15808 


3 D 1 


24489 


25458 


24784 


3 D 2 


24571 


25507 


24802 


3 D 3 


25270 


25591 


24844 


x Pi 


25068 


24990 


25078 



calculated E(M CDF) differ from the experimental data 
in the range 0.3% — 22.0% and the sequence is incorrect 
as 1 P\ lies below the 3 Dj. A correct sequence of the lev- 
els can be obtained by saturating the valence correlation 
effects by including the virtual orbitals 6d and 6/ gener- 
ated as correlation orbitals. The results of the calculation 
are give in column MCDF2 of Tablefv]. The 3 Dj and 1 P 1 
are in good agreement with the experimental data and 
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the 1 P\ lies above 3 Di and 3 r>2 levels. However, it is also 
above 3 D 3 , which is in disagreement with the experimen- 
tal data. This indicates that the correlation orbitals 6d 
and 6/ lower the 3 Dj levels but has little effect on the 
3 Pj and 1 P 1 levels. 

The MCDF calculation captures the important 
valence-valence correlation effects. However, it is not 
suitable to capture the core-valence and core-core cor- 
relation effects, which require a large number of CSFs. 
CI calculations within the CSF space having excitations 
from the core and valence shells to the virtual can capture 
these correlation effects. The virtual orbitals required 
for the CI calculation is generated in layers, where one 
layer is a set of orbitals of s, p, d, f and g symmetries 
having same principal quantum number. Higher angular 
momentum orbitals h and above are not included in the 
calculation. Each layer is generated by an MCDF-EOL 
caculation of the CSFs used in the previous calculation 
and the CSFs obtained by single excitation from the con- 
figurations 6s 2 , 5cffis, 6s6p, 6p 2 , 6p6f and 5d6p. 



TABLE V: The excitation energies of the low-lying levels and 
the magnetic hyperfine structure constants. 





Excitation energies 


Hyperfine 


Constant A 


Level 


Expt 


Theory 


Expt 


Theory 


3 Po 


17288 


15693 

17282" 

17359 6 






3 Pi 


17992 


16370 

17997" 

18089 6 


-1094.0(7) c 


-1065 
-1094" 


3 P 2 


19710 


17989 

19750" 

19836 6 


-738 d 


-747 
-745" 


3 Zh 


24489 


23468 

24441" 

24936 6 


563 c 


601 

596 ° 


3 D 2 


24571 


23567 

24697" 

25180 6 


-362(2) c 


-307 
-351" 


3 D 3 


25270 


23744 

25247" 

25676' 


-430(1) c 


-454 
-420" 


'Pi 


25068 


24430 

25074" 

27271 6 


59 d 


152 
191" 



"See ref 
6 See ref 
c See ref 
d See ref 



m 

f[M] 



The results of a CI calculation within the CSF space 
spanned by all possible excitations from the valence shells 
and single excitations from the core-shells 5s, bp and 4/ 



TABLE VI: 




Fo(bsbp) 


r\ (osop ) 


r\ [osop ) 


1 S (6s 2 ) 




0.51 


4.98 






0.54(8)" 


4.40(80)" 






0.44 6 


4.44(80) 6 






0.549(4) 


a on / o r\N c 

4.89(80) 






0.553(13) / 


4.13(10) / 








4.02 9 








4.26 h 


3 Di(5d6s) 


3.03 


2.64 


0.28 




2.61(10)" 


2.26(10)" 


0.27(10)" 






2.2(1) d 


0.24 b 



"Sec rcf[t 
b See ref [[ 
' Sec ref j 
d See ref[[ 
e Scc rcf[p 
■^See ref[E 
»See ref[E 
h See ref [ 



with respect to the 6s 2 , 5<i6s, 6s6p, 6p 2 , 6p6f and 5d6p 
are given m Table|v|. The results of earlier calculations 
are also given in the table. The sequence of the excita- 
tion energies are correct but the diviation from the ex- 
perimental values is still large. However, the level 1 P\ is 
better than the coupled-cluster results and the hyperfine 
constant of the same level is in better agreement with 
the experimental data. The important reduced matrix 



elements of the electric dipole are given in the Table: VI. 

The £Tpnc calculated using the theoretical results is 
—86.01 x lCr n i(— Qw/N)ea . The important contri- 
butions are from the intermediate states 3 Di(5d6s) and 
1 5 (6s 2 ), each of these yield -90.62 x IQ-^i-Qw /N)ea 
and 24.74 x 10 -11 i(Qw /N)eao respectively. The can- 
cellation due to the 1 S'o(6s 2 ) is large but is nearly 
compensated by the contributions from the interme- 
diate states 3 Si(6s7s) and 3 Pi(6p 2 ), which totals to 
—21.78 x 10~ 11 i(Qw/N)eaQ. However, as mentioned 
earlier, the energy spacing of the important levels are 
not correct. In addition < 3 P (6s6p)||D|| 3 Di(5d6s) >, 
< 3 P (6s6p)|| J D|| 3 S'i(6s7s)> and | < 150||Z?||1P1 > | also 
show variations from the experimental data and previous 
theoretical calculation Q]. Using the experimental exci- 
tation energies, the average of the experimental data of 
| < lS'OHP'lllPl > |, and the other |.D/i|s calculated by 
Porsev et al[|l) gives -12.31 x 10- 10 i(Q w /N)ea . 

The contribution from the intermediate state 
3 P>i(5d6s) is -13.09 x 10- 10 i{Q w /N)ea which 
is calculated using the reduced matrix element of 



3 Po(6s6p)> from the present work, 
| 3 P (6s6p) > from Porsev et alpl] 



< 3 Di(5d6s) 
the < 1 Pi(6s6p)||D 
and the enrgy denominator from the experimental data 



The estimate of the < 3 P>i(5d6s) 



i P (6s6p) > 



using the wave functions we have calculated is expected 
to be better as the magnetic dipole hyperfine structure 



6 



constant of 1 Pi(6s6p) is in better agreement with the 
experimental data. 



that 3 Po(6s6p) — 1 Pi(6s6p) is indeed a promising 
candidate for the -EIpnc measurements. 



IV. CONCLUSIONS 

Our MCDF calculations shows that the 
£iPNc( 3 -Fb(6s6p) - 1 P 1 (6s6p)) PNC transition am- 
plitude of 171 Yb is more than two order of magnitude 
larger than that of the i?lpNc(6s — 7s). Our calculation 
indicates the experimental proposal of Kimball[[l) shows 
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